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The formation of patterns of peaks on the free surface of a ferrofluid subject to a magnetic field 
normal to the undisturbed interface is investigated theoretically. The relative stability of ridge, 
square, and hexagon planforms is studied using a perturbative energy minimization procedure. 
Extending previous studies the finite depth of the fiuid layer is taken into account. Moreover, 
adding the wavenumber modulus k to the set of variational parameters also the wavenumber selection 
problem is addressed. The results are compared with previous investigations and recent experimental 
findings. 
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^\ [ When a ferrofluid layer is subjected to a vertically oriented and uniform magnetic field, above a critical value of 
the field strength a hexagonal pattern of peaks appears on the surface of the liquid. This Rosensweig or normal field 
instability was first observed by Cowley and Rosensweig in 1967 [0j2j. Further increase of the magnetic field up to a 
^ . second threshold gives rise to a transition from the hexagonal to a square planform . 

. ■ The arrangement of peaks resulting from the Rosensweig instability into patterns of different geometry is just one 
S , particular example from an impressive variety of pattern formation in physical systems ||^. It is well known that 
although the instability threshold itself can be obtained from a linearized version of the underlying equations, the 
pattern selection problem requires the inclusion of non-linear terms. A standard procedure to probe the non-linear 
regime perturbatively is by means of amplitude equations 
' However, unlike many other examples of pattern forming physical systems discussed in the literature the surface 
^ profile of a ferrofluid in a static magnetic field is an equilibrium structure. Accordingly the relative stability of plan- 
forms and possible transitions between different patterns can be investigated theoretically by studying the appropriate 
thermodynamical potential. Still the problem is a complicated non-linear task since the local magnetic field deter- 
' mining the surface profile in turn depends on the surface deffection via boundary conditions. As a consequence the 
I variational minimization of the thermodynamic potential in the surface profile cannot be accomplished exactly and 
T-H • one has to resort to approximate procedures applicable in the experimentally relevant situation of slightly overcritical 
magnetic fields. 

In classical investigations along these lines Gailitis |^] and Kuznetsov and Spektor ||] analyzed the stability of 
the different patterns of the Rosensweig instability by means of an energy minimization principle. These as well as 
. related investigations using methods of functional analysis [^,|l^ or a generalized Swift-Hohenberg equation 12 
were confined to fiuid layers of infinite depth. 



On the other hand experimental investigations are usually done with rather thin fluid layers and the effects of 
finite thickness on the linear regime have been studied recently in great detail. For example the dispersion relation of 
^ surface waves on a ferrofluid layer of arbitrary thickness was determined in jl^ and the influence of viscosity on the 
linear dynamics was elucidated in ||l^ . Taking advantage of the dependence of the threshold of instability and the 
wavelength of the most unstable mode on the thickness of the fluid layer in a clever wa y, i t is e.g. possible to measure 
both the normal and the anomalous dispersion branch of surface waves on ferrofluids [|l5| . 

In the present paper we complement these investigations by a thorough study of the weakly non-linear regime of the 
Rosensweig instability slightly above the critical magnetic field for a ferrofiuid of arbitrary depth. Taking the limit 
of infinite layer thickness we also critically discuss the classical findings obtained in |^J^. Moreover we are able to 
quantify the restriction to sufficiently small susceptibilities x of the fluids which was always used in previous studies. 

Our method of investigation is a generalization of the variational minimization of an energy functional already used 
in 1^,11 . Near the instability this functional may be written as a power series in the amplitude of the surface deflection 
and the minimization can be performed explicitly. Moreover our approach allows the theoretical investigation of the 
wavenumber selection problem addressed also in recent experiments Including the wavenumber k into the set of 
variational parameters we determine the dependence of the wavenumber of the patterns on the magnetic field and 
investigate the infiuence of a varying wavenumber on the stability of hexagons and squares. Therefore by a slight 
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extension of our method we are able to study questions which are not easily accessible to other non-linear methods, 
as e.g. amplitude equations. 

The paper is organized as follows. In section th e basic equations are collected and transformed into a form 
suitable for the calculation of the energy. In sectioii |lIl| a perturbation ansatz for the sur face deflection is put forward 
and the important issue of its consistency and region of validity is discussed. Section IV is devoted to the comparison 
of our findings with the classical results of pattern selection in a fiuid with infinite depth. Subsequently in section ^ 
we consider a fcrrofluid layer with arbitrary thickness and analyze the effects of the finite dept h on the stability of 
the patterns. In section VI we address the wavenumber selection problem. Finally section VII contains a summary 
and compares our findings with recent experimental results. 



II. BASIC EQUATIONS 

We consider the situation sketched in Fig. |^. A horizontally unbounded ferrofiuid layer is subjected to an external 
magnetic field Hq which in the absence of any magnetic permeable material is of the form Ho(x,?/, z) = —Hoez- The 
incompressible magnetic fluid of density p, surface tension a and susceptibility x is bounded from below at z = —d 
by an impermeable material and has a free surface described by z = (^{x, y) with the magnetically impermeable air 
above. The gravitational acceleration g = ^ge^ acts parallel to the z-axis. Our aim is to determine which static 
profile C{x,y) develops for a magnetic field Hg strong enough to destabilize the flat surface C{x,y) = 0. 




FIG. 1. Schematic plot of a ferrofiuid layer with infinite horizontal extension in an external magnetic field Ho parallel to the 
gravitational acceleration g. We investigate the pattern formation on the free surface z — ({x,y) for arbitrary depth d of the 
magnetic fluid with susceptibility x > 0- 

Stable configurations of the ferrofiuid surface with infinite horizontal extension are given by minima of the thermo- 
dynamic potential per unit area in the x — y plane 

/[C(^,2/)] = ( fc'(^,2/)- Y J dz HoxH(x, y,z) + a^l + (5.C(x, yjf + {dyCi^, y)? \ ■ (1) 

Here fiQ is the permeability of free space, H(a;, y, z) is the magnetic field in the presence of the ferrofiuid, and the 
brackets denote the average over the x — y plane defined by 

L L 

{Fix,y)) :=^li^^ J dx J dy F{x,y). (2) 

-L -L 

The three terms on the r.h.s. of Eq. (Q) describe the hydrostatic energy, the magnetic energy and the surface 
energy respectively. Note that both the hydrostatic and the surface energy increase when the interface profile starts 
to deviate from the fiat reference state, whereas the magnetic energy decreases. For sufficiently large Hq this gives 
rise to the normal field or Rosensweig instability [|l]j2| . 

The magnetic fields B and H are determined by the static Maxwell equations. 
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together with the boundary conditions, 



V • B = and V x H = 0, 



hm H(a;,?/, z) = -iJoBj, 

z — >±oo 



(3) 
(4) 



[(Ba - Bff ) • n] =0 and [(Ha - Hff) x n] - 0, 



(5) 



[(Bff - Bb) • e,] 



z— — d 



and [(Hff - Hb) x e,] 



z=-d 



(6) 



where n denotes the normal vector on the surface ^(x,y) and the respective subscripts indicate the fields above, in 
and below the ferrofluid. The minimization problem for the thermodynamic potential becomes non-trivial since the 
magnetic field ll{x,y, z) depends on the profile C,{x^y) of the surface. 
Throughout the paper we will assume the linear relation 

B = Aio(l + x)H (7) 

between the magnetic induction and the magnetic field. For the experimentally relevant magnetic fluids and fields 0| 
this approximation deviates about 5% from the exact relation between B and H. 
It is convenient to introduce a scalar magnetic potential "0(2;, y, z) defined by 

H = -VV' (8) 

which by (^) has to satisfy the Laplace equation 

A?A = 0. (9) 
The two characteristic scales of the problem are the critical wavenumber at the onset of the instability. 



for a ferrofluid with infinite depth d — > c» and the corresponding critical magnetic field. 



' (l + x)(2 + x)2Vpg^ 



(10) 



(11) 



which were first derived in Q from a linear stability analysis. Henceforth we therefore measure all lengths in units of 
the capillary length fc^T^, all wavenumbers in units of the critical wavenumber k^ ooi the magnetic field Hp in units of 
the critical value -ffc,oo, the scalar magnetic potential V' in units of -ffc,oo/fcc,oo, and energies per unit area in units of 
a. Moreover it is convenient to introduce the following rescaled magnetic potentials in the space above, in and below 
the ferrofluid respectively: 



(2 + x) 

(l + x)(2 + x) 

HoX 
(2 + x) 



Each of these potentials has to fulfill the Laplace equation (||) . Using the abbreviation 

X 



V 



X 



(12) 

(13) 
(14) 

(15) 



the asymptotic form of the magnetic field for |z| —f oo as specified by Eq. ^ gives rise to the requirements 
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lim dz'4)a^{x,y,z) = - = Van dziphix^y, z) 



2:— »+oo 



(16) 



for the magnetic potentials. Moreover the boundary conditions (|^) and for the magnetic fields translate into the 
following boundary conditions for the potentials: 



dyC, - [dzi^s, - dz^s] 



1+7/ 



0, 



(17) 
(18) 



_ l + r] 



0, 



= 0. 



(19) 



(20) 



z— — d 



Using Eqs. (||) and (|T^) and exploiting the fact that Hq is parallel to the z-axis, we finally get the energy / as 
functional of the surface deflection ( ~ y) in the form 



nc] = ( Y - ^0 {Mz=c - Mz=^d) + ^/^ + (Mr + {dycr) ■ 



(21) 



As stated above, the main problem, rendering a straightforward minimization of / in (^{x, y) impossible, is the rather 
implicit dependence of the potential "^g (x, y, z) on the surface deflection C(a;, y) specified by the boundary conditions 
(0) and (|l|). 



III. THE PERTURBATION ANSATZ 



The variational problem for the energy functional posed in the last paragraph can in general not be solved exactly. 
In order to make analytic progress, we restrict ourselves to the vicinity of the critical magnetic field and assume that 
the amplitude of the surface deflection is still small. It is then possible to expand the energy in this amplitude and 
to retain only the first terms. In this paper we will consider an expansion of the energy up to fourth order in the 
amplitude. The applicability of this approach will be critically discussed below. 

Generalizing the perturbation expansions used in |Q,||,|l9|,^ for the one-dimensional variant of the Rosensweig 
instability, we write the surface profile in the form 
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C(a^, y) = / ^ cos(k„ • r) 



17 



i3„cos(k„ • r). 



(22) 



where r = ix,y) and k — (kx,ky) are two-dimensional vectors. The terms with n = 1,...,4 constitute the main 
modes. Their corresponding wavevectors ki to k4 have all the same modulus k whereas their mutual orientation 
is chosen such that the ansatz allows the description of ridges (e.g. Ai = A, A2 = A3 — A4 — 0), squares (e.g. 
Ai = A4= A, ^2 = ^3 = 0), and hexagons (e.g. Ai = A2 ^ A3 = A, A4 = 0) (see Fig. |). 

The terms with n = 5, 17 are higher harmonics with wave vectors ks to kiy being linear combinations of two 
wave vectors of the main modes and amplitudes _B„ of order 0{A'^). These terms are needed to satisfy the minimum 
conditions for the energy functional to the required order 0{A^). The intuitive meaning of this fact is that reliable 
results on the relative stability of different planforms requires some information on the deviation of the nonlinear 
surface profile from the simple cosine shape describing the linear instability. On the other hand there is no need 
to include even higher harmonics in the ansatz (22) since the corresponding contributions would average to zero in 
the X — J/— integrations in the definition ( pl| ) of the energy functional. In conclusion the chosen ansatz is the only 
consistent perturbation ansatz for the energy which includes terms up to fourth order in the amplitude of the surface 
deflection. The values of the amplitudes Ai, . . . ,A4 and B5, . . . , Bn as well as the wave vector modulus k of the main 
modes are the free parameters which may be used to minimize the energy. In particular the possible minimization in 
fc is a special advantage of our approach since it allows an analysis of the non-linear wave vector selection problem 
not easily accessible to other approaches. 
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FIG. 2. Wave vectors ki to k4 of the four main modes in the perturbation ansatz ( p^ (arrows). k4 is perpendicular to 
ki and the angle between ki and k2 or k2 and ka is 7r/3. The endpoints of the wave vectors of the higher harmonics are 
represented by crosses. Each of these wavevectors is the sum or difference of two of the main wavevectors. 

Because of the boundary conditions at the surface of the fluid the magnetic potential ip must have a similar 
dependence on x and y as C,{x,y). The ansatzes 
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k„|z 
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n=0 



cos(kn • r), 
l''"l^cos(kn -r), 



(23) 



fulfill the Laplace equation ^ and the asymptotic boundary conditions (|T^). Omitting the mode ko (0,0) in the 
potential ■(/'ff fixes the free constant of the potential if). Plugging Eq. ( p3| ) into the four remaining boundary conditions 
Eqs. (p^)-(|2C|) we can determine the parameters u„,w+,ti~ and Wn as functions of the surface deflection amplitudes 
by expanding the equations up to third order in the An. The resulting expressions specify the dependence of the 
magnetic field on the surface profile which is then used to calculate the energy / using Eq. (|T]) up to fourth order in 
the amplitudes A„. The resulting dependence of / on the higher order amplitudes i?„ is simple and the minimization 
in the _B„ can be performed explicitly. Subtracting from the resulting expression the reference value of the energy for 
a flat interface we finally arrive at 



fd{{A.n},k) :-/({A„},fc)-/(C = 0) 



(24) 



Al 



Al]-^{e,k)[AiA2A^] + -g{eMA\ 



A\ 



At 



At\ 



--gn{e,k)[AlAl+AlAl 



AlA\] + -g,(e, k)[AlAl + AlAl] + -.9„(e, k)[A\Al] + 0{A% 



(25) 



This expression gives the energy of a ferrofiuid layer of depth d with surface profile C,{x,y) as specified by (|2^) at 
arbitrary strength of the external field Hq up to fourth order in the defiection amplitudes A„. The dependence on 
Hq is expressed via the supercriticality parameter 



Hi 



- 1. 



(26) 



The function Z(e, k) is given by 



l{e,k) 



1 



-2kd 1 

^k-^e. 



(27) 



It fixes the critical wavenumber fcc,d for the onset of the instability and the corresponding critical magnetic field iJc.d 
for a ferrofiuid with finite depth d via the conditions 
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^(^ = 0)L=fe„, =0 (28a) 



and 

dlje = 0) I 
dk lfc=fcc 



0. (28b) 



Determining the minima of / with respect to the amplitudes An, one finds that e — > with A„ ^ as expected: the 
external field will be near to its critical value if the equilibrium values of the amplitudes An are small. For ridges and 
squares the explicit scaling is e ^ A^ for A — > 0, in the case of hexagons one finds e ^ —A for small A. This in turn 
implies that to the desired order O(A^) of the expansion of / in the amplitudes A„ only the e-dependence of I and 7 
has to be retained. The exact form of these dependencies is l{e, k) — Z(0, k) + edl{e, k)/de and 7(6, k) = {1 + e)^{k). 
In this way we arrive at our main result for the energy /(i({A„}, A;) of the surface defiection which describes the 
equilibrium configurations consistently up to fourth order in the amplitudes An'- 

U{An}, k) = -hie, k)[Al + AI + AI + Al] - (1 + e)^{k)[A,A^A^] + \g{k)[A\ + A\ + Al + Af, + 



+ \gn{k)[A\Al + AlAl + AlAl] + \gt{k)[AlAl + AlAl] + \gn{k)[A\Al] + 0{A% (29) 



The explicit expressions for the various coefficients in (|29|) are rather long. We therefore display their form here only 
for the comparatively simple situation of a fluid of infinite depth, d — > (X): 

7oo(fc) = |r,fc2, (30^) 
..(^) = ^^^(8-3.)-^-|!!^, (30b) 
,.,.(fc) ^ h^HW - 7.^V3 + 3V3 - - 3) - ^ f_^',%'^^^ (30c) 



,o,{k) = fc3(^V3\/2 - 1 + 47?2 - lifViV^- ^k) + 
4 4 16 

1 T/^fc4(18 + \/3 - 4V2 - 4\/3%/2) 1 7?^fc4(18 - V3 + Ay/2 - 4\/3V2) 
8 2 - fc\/3\/2 + fc2 + A:V2 - VS 8 2 - A:V3\/2 + F - fcV2 + VS ' 



(30d) 



IV. FERROFLUID LAYER WITH INFINITE DEPTH 

In this section we consider a ferrofiuid of infinite depth d 00 and fix the wavenumber modulus k at its critical value 
^c.oo = 1- This enables us to compare our findings with the well-known results of previous theoretical investigations. 

Possible surface patterns together with their amplitudes are given by stationary points of the energy function 
foo{{An},k = 1). We will in particular study ridges ("rolls") given by Ai — Aji with A2 = A3 = A4 = 0, squares 
represented by Ai = A4 = As with A2 = A3 = 0, and hexagons described by Ai = A2 = A3 = Ah with A4 — 0. The 
last planform is shown in Fig. Using Eqs. (^^, (|2^) and (^0|) we find for the corresponding amplitudes 

(31) 




Me) = J ■ , (32) 
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and 



, , , 1 7oo(l + e) + V7i(l + + ^e{2gh 

AH[e) = 

^ 'igh,oo + goo 



(33) 



respectively. Together with x slso 7 is always positive]^. Hence only up-hexagons with a peak instead of a dip in the 
center are stable. 

Our results for the hexagons differ slightly from the classical ones reported in |Q in which the dependence of the 
cubic term in Eq. ( p9[ ) on e is missing. In |^,^ this dependence was omitted because the discussion was restricted to 
fluids with 7 <C 1 so that the approximate scaling e ^ is very similar to the true one e ^ —jA + O(A^). Otherwise 
our results coincide with those of Gailitis, in particular using k — fee. 00 = 1 in ( pO| ) the respective expressions found 
in 1^ are reproduced. The investigations of Kuznetzov and Spektor |^ were restricted to first order in rj and are 
contained in both our results and those of Gailitis. 
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FIG. 3. Lines separating regions in the X ~ ^ plane in which the amplitudes of ridges Ar, squares As and hexagons Ah 
are smaller than 0.25 (full lines) or smaller than 1 (dashed lines). The thin dotted line represents the results given in Hi. For 
susceptibilities larger than the critical values XR^XS ^^'^ Xh the amplitudes diverge in the considered approximation (E2|). 

Before addressing the relative stability of the different planforms introduced above it is important to characterize 
the domain of validity of our central expression ( |29| ) for the energy. Being a perturbative result the corresponding 
amplitudes An must not be too large. Fig. |^ gives a quantitative impression of what that means by displaying lines in 
the x~ ^ plane corresponding to fixed values of the amplitudes of different patterns. It is in particular important to 
realize that for sufficiently large values of the susceptibility x the fourth order coefficients in the energy functional may 
change sign. The corresponding patterns than appear through a backward bifurcation |p^ . At the same time higher 
order terms (i.e. sixth or higher order in the amplitudes A„) in the energy are necessary to saturate the instability. 
The concrete values of x for the planforms considered are XR ~ 2.54 {goo becomes negative), xs ~ 0.41 (gn.oo + .9oo 
becomes negative), and xh ~ 1-05 {2gh^oo +5oo becomes negative). Beyond these values of x our ansatz (p2) is unable 
to describe the arising pattern. The value of XR fcst calculated in ||l9[] . 

Keeping this limitation of our approach in mind we have investigated the pattern selection problem by study- 
ing the character of the extremum of the energy functional (|29| ) at the ridge, square and hexagon solutions 
given in Eqs. (pl|)-(|3^. This is rather similar to a linear stability analysis of the fixpoint solutions of the cor- 
responding amplitude equations The results are as follows: Ridges are never stable, since gn,oc < 9oo al- 
ways. We have always g^o ~'r 9n,oo < 5/1, 00 + gt,oo and g„,oo < 5oo- Hence (for x < Xs) squares are stable if 



^This holds true also in the case of a finite layer depth d. 
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(,goo + gn.oa - gh,oc - 5t,oo + V(5oo + 5n,oo - 5/>,oo ^ gt.ooY ^ ^llc{9<=° + 9n,oo)T t J f 

e > es — — „ , r . in order tor hexagons to 

474(500 +5«,oo) 

be stable (for x < Xh) the foUowing three conditions have to be fulfilled. 

,ooH" ,900 )(2.9. 



• e > Ei? 



7oo ,00 + 5oo)(2gh 



• ^ 7^ .9h,oo + 2goo 

• either e < = — ^ 77 r or 5,1.00 < 5oo, 

(74 + 9h.oo - 9oo){gh..oo - 9oo) 

• either e < eh = . „ ^ ^ ■ rr^i ■ r Or 25,,.oo + ffoo < 2gt.oo + 5n,oo- 

(7^+ ^gh.,oo+ goo- 25t,oo- 5n,oo)(2g/i,oo+ ffoo- 2gt,oo- ffn,ooj 

These results are displayed in Fig. ^. For regions in which two different patterns are stable, i.e. in which different 
local minima of the energy exist, we have also determined the respective Maxwell values of e at which these minima 
have the same value. For e = ens the energy of the hexagon pattern is equal to the energy of the square pattern, for 
e = ej;^^' it is equal to the energy of the flat surface. 

For all susceptibilities either the relation gh,oo < 5oo or the relation 2gh 00 + 5oo < '2gt,oc + ffn.oo is violated and 
consequently hexagons always become unstable for sufficiently large e. 

Of particular importance is the line denoted by eh- For e > eh the hexagon pattern is unstable to squares which in 
turn for x > Xs £^re not saturated by the fourth order term in the energy. In that case within our perturbation ansatz 
we are not able to predict which pattern will show up. For x > Xo ~ 0-56 this transition occurs right at onset and 
our analysis of hexagons is therefore restricted to % < 0.56. This stability criterion for hexagons was to our knowledge 
not discussed in the literature before. It was not found by Silber and Knobloch since with their approach the 
relative stability between hexagons and squares could not be determined. It was also missed by Gailitis Q and by 
Kuznetsov and Spektor M since they considered only the case x ^ 1- 




FIG. 4. Limits of the stability regions (solid lines) and Maxwell points (dashed lines) in the x~ ^ plane. The left figure refers 
to positive values of e, the right to the subcritical region e < 0. Squares are stable if e > eg. Hexagons are stable if e < en, 
t < th and € > ep. The two thin dotted lines illustrate exemplarily the corresponding findings if the dependence of the cubic 
term in Eq. (^) on t is neglected. 



V. FERROFLUID LAYER WITH ARBITRARY THICKNESS 



The aim of this section is to investigate how the finite depth d of a ferrofiuid layer affects the pattern selection. 
Only ferrofluids with relatively small susceptibilities x < Xs will be considered in order not to leave the region of 
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applicability of our perturbation approach (see the discussion above). 

From the linear stability analysis |l|,|2l[ it is well known that a thin layer retards the onset of the instability and 
shifts the wavenumber of the unstable mode to larger values. The dependence of the critical magnetic field -ffc,d and 
of the corresponding critical wavenumber fcc,d on the layer thickness d can be determined from Eqs. (^8|). It is shown 
in Fig. 1^. If the depth d of the fluid is larger then the critical wavelength at infinite depth Ac,oo — Svr/fcc.oo, which is 
typically of the order of 10 mm, the finite thickness of the layer can be ignored. Even if the layer is very thin, the 
parameters characterizing the linear instability of the flat surface are modified only within a few percent for fluids 
with small susceptibilities. 




FIG. 5. Critical magnetic field He and corresponding critical wavenumber kc as functions of the layer thickness d in units of 
Ac,tx3, where Ac,oo is the critical wavelength for the layer with infinite thickness. Both critical values are altered less then 2% 
for a fluid with x = 0.05 (solid lines) and less then 7% for a fluid with x ~ 0.30 (dashed lines). 



To start with the non- linear analysis we have first to determine the dependence on d of the critical values XRtXS 
and xh beyond which higher order terms in the energy functional are necessary. The corresponding results are shown 
in Fig. H All values of the critical susceptibilities get eventually smaller if the layer depth decreases. They do not cross 
each other such that xs is the smallest one for all values of d. We have included the result for ridges also although 
these are unstable in the setup of Fig. [l] for all d and x- However, ridges can be stabilized by an additional magnetic 
field tangential to the undisturbed surface [22[ and therefore our results give first informations on the behavior in this 
modified experimental setup. 

For very shallow layers the amplitude of the arising pattern may become larger than the layer depth already at 
onset. As a consequence the ferrofiuid layer disintegrates into unconnected regions and "dry spots" occur. This 
topological change is outside our theoretical analysis. The corresponding value of x denoted by Xd is also shown in 
Fig. ||. As expected it becomes the smallest of all critical x- values for d ^ 0. 
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FIG. 6. The critical susceptibilities xn, xs, and xh beyond which higher order terms are necessary in the energy functional 
in order to describe the respective pattern consistently as function of the layer thickness d. Also included are the values for 
XO beyond which hexagons are unstable to subcritical squares introduced in section ^ and x > Xd{d) for which the pattern 
amplitude at onset already exceeds the layer thickness. In both cases our theoretical model ceases to be applicable. 



In order to investigate the stability of the patterns arising on a layer with finite thickness we can perform basically the 
same analysis as in the previous section. We just have to fix the wavenumber k to its critical value /cc^d corresponding 
to the given thickness d and to use the appropriate value for Hca hi the function k — fcc,d) (see Eq. (p7|)). Moreover 
we have to take into account the dependence of the coefficients 7, gh, gt and gn on the depth d. In this way we can 
finally determine for any given depth d the stabilityregions of the different patterns by means of the corresponding 
energy function f{{An}, k — fcc,d) as given by Eq. (|29|). 
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FIG. 7. Limits of the stability regions (solid lines) and Maxwellpoints (dashed lines) as functions of the layer thickness d for 
a ferrofluid with x ~ 0.05. Ridges are always unstable, squares are stable for e > eg and hexagons are stable for ej? < e < en- 
In the gray region the instability results in the formation of dry spots right at onset, so that our theory no longer applies. The 
lower part of the figure describes the subcritical region. 



10 



The results for the various stabihty boundaries and corresponding Maxwell points for a fluid with susceptibility 
X = 0.05 are shown in Fig. |^. All the critical values of e at which transitions between different planforms arise 
increase in absolute value with decreasing layer thickness d. Simultaneously the hysteretic behaviour becomes more 
pronounced. As in the linear analysis we find that the effects of finite depth can be neglected, as long as the fiuid 
layer is thicker then the critical wavelength Ac,oo- However, for a thin layer the stability of the hexagons and squares 
is changed noticeably. For instance the value of ch at which the transition from hexagons to squares takes place may 
increase by a factor of 3 implying an increase of the magnetic field of about 10%. As discussed above for very thin 
layers our theory is no longer applicable since the dips of the hexagonal pattern at onset already reach the bottom of 
the fluid. The corresponding values of d are indicated by the shading in Fig. M. 




-0.1 I —J ^ 
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FIG. 8. Same as Fig. ^ for a ferrofluid with x = 0.30. With decreasing depth d of the layer both xs{d) and xo[d) become 
smaller than the susceptibility of the fluid higher order terms in the energy would be needed to determine the pattern stability. 

Fig. ^ displays the analogous results as obtained for a ferrofluid with the larger susceptibility of x = 0.30. Whereas 
the situation for e < is rather similar to the one shown in Fig. |^, the behaviour for positive e is qualitatively different. 
This is due to the fact that although xs{d — > oo) > x with decreasing depth d both xs{d) and xo{d) eventually 
become smaller than 0.30 (cf. Fig. This gives rise to the divergencies in the energy discussed above and changes 
the stability chart accordingly. 



VI. WAVENUMBER SELECTION 



In many pattern forming systems the wavelength of the first unstable mode is often a good estimate for the typical 
length scale of the emerging pattern. If a whole band of modes is unstable the one with maximal growth rate is 
likely to dominate the arising structure Q. On the other hand the developed pattern is largely determined by the 
non-linearity in the problem and therefore its length scale may also substantially deviate from the relevant scales of 
the linearized theory. 

A particularly gratifying feature of our variational approach to pattern formation in the normal field instability 
in ferrofluids is the possibility to include the wave number modulus k into the set of parameters varied to minimize 
the energy functional. This allows to directly determine the dependence of the pattern periodicity on the external 
magnetic field and to investigate the influence of the variation in wavenumber on the stability of the different patterns. 

We therefore investigate in the present section the wavenumber selection problem for the Rosensweig instability. For 
simplicity we will consider a ferroffuid with infinite depth only. As shown in the previous section this is a very good 
approximation as long as d > Ac,c!o- We will compare our findings with recent experiments in which this condition is 
met. 



11 



The necessary calculations are similar to those reported in section IV with the only extension that we no longer fix 
the wavenumber modulus k in advance but include it into the set of parameters with respect to which the energy is 
to be minimized. 
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FIG. 9. Surface deflection ("(0, 0) at the cusps of the pattern and emerging wavenumber k in dependence on the supercriticality 
parameter e. For a fluid with x ~ 0-35 squares are only stable in the region eg < e < ea. Ridges are always unstable. 

A typical result is shown in Fig. ^displaying the maximal surface deflection (,{0, 0) and the corresponding wavenum- 
ber k of the patterns. As can be seen from the figure, the wavenumber of the patterns decreases with increasing 
supercriticality. This is in particular pronounced for the square pattern and implies that in this case the system may 
lower its energy very efficiently by increasing the wavelength of the pattern. 

This effect introduces an additional threshold for e into the stability chart shown in in Fig. If e > eA we 
find k and the pattern disappears by increasing its wavelength to infinity. This process is accompanied by an 
unbounded increase of the amplitudes Ai , A4 which carries us outside of the validity of our perturbation ansatz ( |2^ ) . 
The stability regions of the other pattern remain qualitatively the same as without minimization in k as can be seen 
from a comparison between figures ^ and |^. 
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FIG. 10. Limits of the stability regions (solid lines) and Maxwell points (dashed lines) in the X~ ^ plane if variations of the 
wavenumber modulus k are taken into account. Squares are stable if es < e < ex. Hexagons are stable if < e < min(e£f , e^). 
For large values of x the stability region for squares gets smaller due to the wavelength instability occurring for e > e\. 
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For supercritical magnetic fields a whole band of wavenumbers giving rise to stable patterns exists. For hexagons 
this stability region is shown in Fig. |ll|. At fixed e the hexagonal planform becomes unstable if the wavenumber 
is increased sufficiently. In |^ such an increase of the wavenumber was achieved experimentally by compressing the 
hexagonal pattern in a hopper. The observed transition from a hexagonal to a square pattern with smaller wavenumber 
occurs when the limit Ch of the stable fc-band is exceeded and is hence in qualitative agreement with our theory. 




FIG. 11. Stability band of tlie hexagonal pattern in the k — e plane for a ferrofluid with x = 0.35. Hexagons become unstable 
when e > eh- The marginal stability curve (dashed line) is also displayed. 

VII. DISCUSSION 

In the present paper we have theoretically investigated the formation of patterns on the free surface of a ferrofluid 
resulting from the Rosensweig instability. For a ferrofluid layer of arbitrary depth d in an external magnetic field Hq 
the equilibrium surface profile was determined by minimizing the appropriate thermodynamic potential. We have 
assumed a linear dependence of the magnetization on the magnetic field which is at most 5% off the correct values 
for the experimentally situations with which we compare our findings. The analysis is restricted to the vicinity of 
the critical magnetic field at which the flat interface becomes unstable since it uses a perturbation expansion for the 
energy up to fourth order in the amplitude of the surface deflection. Our findings generalize the classical results of 
Gailitis and Kuznetsov and Spektor ||^ to layers of finite depth d < oo. 

We have found that, at least for the slightly supercritical magnetic fields considered, the finite thickness of the 
ferrofluid layer can be neglected as long as it remains larger than the critical wavelength of the linear instability which 
in typical experiments is about 1 cm. For thinner layers nonlinear effects become increasingly important. This shows 
up quite generally in more pronounced hysteresis effects for the transition between different planforms. We also find 
smaller values of the susceptibilities for which squares and ridges may appear through backward bifurcations and 
correspondingly higher order terms in the energy are necessary to saturate the instabilities. The critical values of 
the susceptibility x beyond which our theory breaks down are displayed in Fig. ^. Together with the findings shown 
in Fig. ^ they quantify the qualitative statement % 1 used in [0,1). An extension of our analysis including higher 
orders similar to what has been done in for the one-dimensional situation seems feasible but rather tedious. 

It is possible to understand the enhanced hysteretic behaviour for thin layers qualitatively. A thin layer inhibits 
surface deformations resulting from the Rosensweig instability (see Fig. ||) since the magnetic field energy is suppressed 
due to the finite depth of the layer. On the other hand for the developed pattern the magnetic fiux is concentrated at 
the peaks of the surface so that the effective layer depth becomes larger. Correspondingly the pattern is additionally 
stabilized and decays only after a reduction of the field stronger than expected from the linear theory. A similar 
reasoning explains why backward bifurcations of ridges and squares are facilitated in thin layers. 

In addition to determining the stable planforms by minimizing the energy functional in the amplitudes of the 
different modes present in our ansatz for the surface profile, we have also addressed the problem of wavenumber 
selection by including the wavenumber modulus k into the set of variational parameters. We found quite generally 
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a decrease of the wavenumber of the pattern with increasing magnetic field. This is in quahtative agreement with 
recent experiments by Abou et al. Q in which the the field intensity was abruptly increased to overcritical values 
and the corresponding wavelength A = 27r/fc of the resulting hexagonal pattern was found to increase with increasing 
field. Likewise it was found that the wavenumber k of the hexagonal pattern is already for slightly supercritical 
magnetic fields smaller then the critical wavenumber kc. The corresponding difference Afc is also found in the theory 
and shown in Fig. ^. Unfortunately a quantitative comparison with the experimental results is impossible since in the 
experiment a ferrofluid with a susceptibility of 1.4 > xh was used. For such a susceptibility the theoretical analysis 
would require higher order terms in the perturbation expansion of the energy. Nevertheless we found that Afc grows 
with increasing susceptibility such that the experimentally observed value of a few percent agrees well with our result 
for Afc at smaller susceptibilities x < Xh- 

Moreover the theoretical analysis has shown that when including the wavenumber into the set of variational pa- 
rameters the square pattern gets rather susceptible to an unbounded increase in the wavenumber. This gives rise to 
a reduced stability domain for squares and may be related to the fact that in the experiments using jumps in the 
field intensity always hexagonal arrays of peaks were found [Q. Quite generally our results verify that the resulting 
planform may depend on the details of the experimentally chosen way to reach overcritical magnetic fields. 

Let us finally stress that the variation of the wavenumber of the patterns with the external field as found here 
theoretically and observed in the experiments does not invalidate our ansatz with a fixed wavenumber used before. 
As reported in Q the wavenumber of the arising hexagonal pattern coincides with fcc if the magnetic field is increased 
and decreased in a quasi-static way. This is probably due to the fact that the boundary conditions suppress the 
emergence or disappearance of peaks in a developed pattern. Therefore using the quasi-static process a metastable 
pattern with kc is produced whereas a jump in field intensity gives rise to the most stable pattern corresponding to a 
smaller wavenumber. 
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